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ABSTRACT 

Aims. We investigate photon quenching in compact non-thermal sources. This involves photon-photon annihilation and 
lepton synchrotron radiation in a network that can become non-linear. As a result the 7— ray luminosity of a source 
cannot exceed a critical limit that depends only on the radius of the source and on the magnetic field. 
Methods. We perform analytic and numerical calculations that verify previous results and extend them so that the basic 
properties of photon quenching are investigated. 

Results. We apply the above to the 2006 TeV observations of quasar 3C279 and obtain the parameter space of allowed 
values for the radius of the emitting source, its magnetic field strength and the Doppler factor of the flow. We argue 
that the TeV observations favour either a modest Doppler factor and a low magnetic field or a high Doppler factor and 
a high magnetic field. 

Key words, gamma-rays: general - acceleration of particles - radiation mechanisms: non-thermal 
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1. Introduction 

Astrophysical 7— ray sources are often subject to the so- 
called 'compactness' problem: when the flux at longer wave- 
lengths is combined with the dimensions of the source, 
which are inferred from light-crossing time arguments, a 
high optical thickness to photon-photon annihilation and 
effective absorption of 7— rays are implied. 

The possibility that high-energy photons could pair- 
produce on soft target photons instead of es c ape i n 
compact sources was first discussed by iJellevI ()1966ll . 
while the ratio of the luminosity to the size of the source, 
L/R, emerges as a determining factor of whet her or not 
a hig h-energy photon will actually be absorbed (jHerterichl 
11974 ). This was followed by work that took into account 
photon-photon annihilation not only as a sink of 7— rays, 
but also as a source of elect ron-positron pairs inside 
non-the rmal c ompact sources fBonometto & Rees' Il97l|; 



Guilbert et al.l 1983; Kazanas 1984: Zdziarski & Lightman 



19851: ISvenssod 1198711 . The aim of these models was 
. to calculate self-consistently the photon flux escaping 
from the sources by taking into account the energy 
redistribution caused by the 7— annihilation-induced pair 
cascades. The assumption was that high-energy particles 
or 7— rays were injected uniformly in a source that also 
contains soft photons, and the system was followed to 
its final steady-state through the solution of a set of 
kinetic equations describing the physical processes at 
work. Several groups have also developed numerical 
codes for the computation of time-dependent solutions 
to the kinetic equati ons taking p hoton-phot on annihi- 
lation into accoun t (ICoppil 19921: [M astich iadis fc Kirfl 
119951: IStern et al.l 119951: lB6ttcher fc Chiand l2002ll . 



These algorithms are commonly used in source mod- 



elling dMastichiadis fc KirkI 119971: iKataoka et al.l 120001: 



iKonopelko et al.ll2003t iKatarzvhski et al.ll2005i r 

One of the assumptions of these models, as stated above, 
is the presence of soft photons in the source, which serve as 
targets for the 7— ray annihi lation. A differen t appr oach has 
been recently presented by IStawarz" fc KirkI (|2007il . hence- 
forth SK. These authors focused on the non-linear effects in- 
duced by photon-photon annihilation and investigated the 
necessary conditions under which 7-ray pho t ons c an cause 
runaway pair productioi£]. IStawarz fc KirkI ()2007D showed 
that there is a limit to the 7— ray luminosity escaping from 
a source, which does not rely on the existing soft photon 
population but is instead a theoretical limit depending only 
on parameters such as the source size and its magnetic field 
strength. Violation of this limit leads to automatic quench- 
ing of the 7— rays. This involves a network of processes, 
namely photon-photon annihilation and lepton synchrotron 
radiation, which can become non-linear once certain crite- 
ria are satisfied. In this case electron-positron pairs grow 
spontaneously in the system and the 'excessive' 7— rays are 
absorbed on the synchrotron photons emitted by the pairs. 
As a result the system reaches a final steady state where 
the 7— rays, soft photons, and electron-positron pairs have 
all reached equilibrium. Therefore this can occur even in 
the hypothetical case when there are no soft photons in the 
source, at least initially. 

Indeed, photon quenching has some interesting implica- 
tions for emission models of 7— rays because it gives a ro- 
bust upper limit on the 7— ray luminosity. The aim of the 
present paper is to investigate the implications of this net- 
work for the high - energy compact astrophysical sources. 
In §2 we will expand the analytical approach of SK, who 
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^ We note that in the same framework an analogous study 
where ultrarelativistic proton s are the primary injected pa rticles 
has earlier been presented bv lKirk fc Mastichiadis! (| 19921 ). 
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studied the dynamical system of photons and relativistic 
pairs using (5— functions for the photon-photon annihilation 
cross section and for the synchrotron emissivity, while they 
treated synchrotron losses as catastrophic. In §3 we will 
use a numerical approach that will allow us to study the 
properties of quenching using the full cross section for 77 
absorption and the full synchrotron emissivity. As an ex- 
ample, in §4 the above will be applied to the 2006 MAGIC 
TeV observations of quasar 3C279 to extract a parameter 
space of allowed values for the source parameters, i.e, the 
radius R and the magnetic field strength B, as well as for 
the Doppler factor 5 of the flow. Finally in §5 we conclude 
and give a brief discussion of the main points of the present 
work. 



2. Analytical approach 

2.1. First principles 

We begin by assuming a spherical source of radius R con- 
taining a tangled magnetic field of strength B. We also as- 
sume that monoenergetic 7— rays of energy (in units of 
rrieC^) are uniformly produced by some unspecified mech- 
anism throughout the volume of the source. If these are 
injected with a luminosity L™^ , one can define the injected 
7— ray compactness as 



£inj _ 7 



(1) 



where ctt is the Thomson cross section. Without any 
substantial soft photon population inside the source, the 
7— rays will escape without any attenuation in one crossing 
time. However, as SK showed, the injected 7— ray compact- 
ness cannot become arbitrarily high because if a critical 
value is reached, the following loop starts operating 

1. Gamma-rays pair-produce on soft photons, which can 
be arbitrarily low inside the source. 

2. The produced electron-positron pairs cool by emit- 
ting synchrotron photons, thus acting as a source of soft 
photons. 

3. The soft photons serve as targets for more 77 inter- 
actions. 

There are two conditions that should be satisfied si- 
multaneously for this network to occur: The first, which is 
a feedback condition, requires that the synchrotron pho- 
tons emitted from the pairs have sufficient energy to pair- 
produce on the 7— rays. By making suitable simplifying as- 
sumptions, one can derive an analytic relation for it - see 
also SK. Thus, combining (1) the threshold condition for 
77-absorption e-yEo — 2, (2) the fact that there is equipar- 
tition of energy among the created electron-positron pairs 
7p = 7e = 7 = 6^,72 and (3) the assumption that the re- 
quired soft photons of energy eq are the synchrotron pho- 
tons that the electrons/positrons radiate, i.e., eo — 67^ 
where b = B/B^^it and Bcrit = {mlc^)/{en) ~ 4.4 x 10" G 
is the critical value of the magnetic field, one derives the 
minimum value of the magnetic field required for quenching 
to become relevant 



B,= 



i^^'^Bcrit- 



(2) 



Thus for B > Bq the feedback criterion is satisfied. Note 
that this is a condition that only contains the emitted 
7— ray energy and the magnetic field strength; moreover. 



it can easily be satisfied, at least if one is to us e the val- 
ues inferred from typical modelling of the sources (jBottcheil 
2007; B5ttcher et al. 2009). 

Along with the feedback criterion, a second criterion 
must be imposed for the quenching to be fully operative. A 
simple way to see this is with the following consideration. 
Assume that the 7— rays pair-produce on some soft photon 
and that the created electron-positron pairs cool by emit- 
ting synchrotron photons. Because an electron emits several 
such photons before cooling the critical condition occurs if 
the number density of the 7— rays is such that at least one 
of the synchrotron photons pair-produces on a 7— ray in- 
stead of escaping from the source. Thus the condition for 
criticality can be written as 



Afs{e-f /2)n{ej)a~^~f R > 1, 



(3) 



where n{e-y) is the number density of 7— rays, a^j is the 
cross section for 77 interactions and Afs{'l) is the number 
of synchrotron photons emitted by an electron with Lorentz 
factor 7 before it cools. Assuming that the pair-producing 
collisions occur close to threshold, approximating the cross 
section there by a-yj ~ (7^/3 and using Msi'y) — and 

K' R 



the critical condition can be written as 



(4) 

where the feedback condition ([2]) and eq. ([T]) were also 
used. Taken as equalities, relations ^ and (U) define the 
marginal stability criterion, which essentially is a condition 
for the 7— ray luminosity. 

2.2. Kinetic equations 

We begin by writing the kinetic equations that describe the 
distributions of 7-ray photons, soft photons, and electrons 
in the source. These are respectively 



dr 



dno{x,T) 
dr 
and 

dr 



no = Ql 



^11 



= q: 



■77 



c. 



syn I 



(5) 



(6) 



(7) 



where n, np and, rig are the differential 7-ray, soft photon, 
and electron number densities, respectively, and e^, x, and 
7 are the corresponding energies normalized in rrioC^ units. 
The densities refer to the number of particles contained 
in a volume element a^R. In other words, if hi expresses 
the number of particles of species i per ergs per cm^ , then 
ni — hi{aTR){meC^). Time is normalized with respect to 
the photon crossing/escape time from the source, ta- = R/c. 
Thus, r, which appears in the kinetic equations, is dimen- 
sionless and equals t = ^. The operators Q and £ denote 
injection and losses, respectively. The only processes that 
we take into account are 77 annihilation and synchrotron 
cooling of the produced pairs. The synchrotron emissivity 
is approximated by a 5-function, i.e., js{x) = jo6{x — eo), 
where eg = ^7^ is the synchrotron critical energy. In all 
cases we will assume a monoenergetic 7— ray injection at 
energy e^. 
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2.2.1. 6 - function approximation for the a-f^ cross section 

First we will examine the stability of the system using the 
simplest form for the photon-photon annihilation cross sec- 
tion (Zdziarski & Lightman 1985.) : 



(x) 



xd 



(8) 



where x, try are the energies of the soft and 7-ray photons, 
respectively. 

The particle injection and loss operators take the fol- 
lowing forms: 



ri 



and 



^l^h-^''x-''^n,{^h^T) 

3 6'-y 



(9) 
(10) 
(11) 
(12) 

(13) 



because each photon-photon annihilation results in a pair 
of leptons, each one with approximately half of the initial 
7-ray energy. In the above equations the 'magnetic com- 
pactness' was introduced 



(14) 



The system of eq. ©-Q can now be written as 
dn{e-y) 



dr 

dno{x) 
dr 

dr 



_^^^^^_2no(2/e,)n(e,)^^^^ 

3 6'-Y 



^no(x) 



27) 



4 no {2/€^)n{ej) 



t^^^ (7'"e(7)) 



<5(7-e^/2)-f 



(15) 



In this section we will examine the stability of the trivial 
stationary solution of the system 



(16) 



which corresponds to the free propagation of 7-ray pho- 
tons through the source. To investigate the stability of the 
system, we assume that initially arbitrarily small pertur- 
bations of the soft photon Uq and electron n'^ densities are 
present, which lead to the perturbation of the 7-ray photon 
density n' . After linearization, the system (jlSp becomes 



dn' 

dr 
dK 
dr 



27) 



X i/2nl 



4 rinn 



+ ^^^(7-^7/2) + |/B^(7X) 



(17) 



The stability analysis of the system above can be simpli- 
fied when one works with the Laplace transformed number 
densities: 



ng(x,T) no(x,s) 



dr ^0(2;, r)e 



with s being the solution to the eigenvalue problem: 
(. + l)n'(.„.) = -^<L(^/±^ 



{s + l)n'o{x,s) = -lBb-^^^x-^/^n'^{^/^,s) 



s<(7,s) = +- ' ^ J(7-e^/2) + 
+ ^/b^ (7'<(7,s)) • 



(18) 



(19) 



(20) 



(21) 



For large r both the soft photon and electron distributions 
behave as e'*'^. If s = — 1, then Uq = n'^ = and the 7-ray 
photons freely escape from the source. If s > 0, one finds 
solutions that grow with elapsing time. The marginally sta- 
ble solution, which is obtained when one sets s = 0, will be 
examined below. 

Combining equations (Tinil - d?!]) leads to an ordinary dif- 
ferential equation for the perturbed electron distribution 
with the following solution: 



7^<(7)-C- 



where jci 
one finds 



2?^<(7c 

3fe3/2,/2^ 



^. Setting 7 



<(7) 



7 



7ci + 



2n 



(22) 

7cr in the above equation 
(23) 



363/2^/21 



Then the perturbed photon distribution is found using 
eq.®- 

As already mentioned in §2.1, the 7-ray compactness 
cannot become arbitrarily high but instead saturates at a 
critical value When << only a negligible frac- 
tion of the hard photon luminosity will be absorbed, i.e.. 



— >■ 0, or equivalently oc Z™-". We call critical the in- 
jected luminosity required for l!^ — to occur. After this 
point any further increase of the injected luminosity will be 
seen as an increase of the luminosity in soft photons. Using 
eq.dlll), (gni, and dM]) one derives 



t_,/t V'^7' S,ci 



where Xn 



spect to we find 



dx xn'g{x;C'^^), 



(24) 



bei/A. Solving the above equation with re- 



7,cr 



(25) 



To compare our results with the value l^^ ~ 3 found 
by SK, who have used catastrophic synchrotron losses, we 
have to take the following steps: 



4 



Petropoulou & Mastichiadis: Automatic 7-ray quenching in compact sources 



■3.2- 



■3.4- 



M^3.6- 




■3.8 



■4.2- 



ton losses becomes 
= "(27) 



B=40G 
B= 4 G 



(27) 



Following the same steps as those described in §2.2.1, we 
find for s = 



where A — and 7max = £7/2- Equation 



(28) 



can be 



solved iteratively when writing n'^{'^) as a sum of approx- 
imations, i.e., ng(7) = n*^"^ + rS'^'^ + • • •. Assuming for the 
first approximation the form n'^^ — C7~^, we find the other 
terms of the series expansion: 



+' 



,(2) 



2 \2 



4.2 



4.4 



4.6 



4.8 



5.2 



5.4 



5.6 



(29) 
(30) 



Fig. 1. Log-log plot of the critical luminosity compactness where A = ^ {-^ 4-j and |A| < 1 for typical values 



as a function of the injected 7- ray photon energy e^, 
wiien the 5-function approximation for the cross section 
a^~f was used for two values of the magnetic field strength: 
B = 40 G (solid line) and B = 4 G (dashed line). The 
values calculated by SK using only catastrophic losses are 
marked with crosses. 



1. The case of synchrotron cooling degenerates to the 
'catastrophic losses case' when the threshold energy 
for absorption 2/e^ equals the maximum soft photon 
energy a;,„ax- 

2. In SK the volume element is defined as V = ttR^, 

whereas we used V — 4/37ri?'^. Thus, the replacement 
l^^ I^SK gi^o^ij ^^^^ 

3. Finally, in SK the luminosity compactness is not given 
by the conventional definition, which we adopted, but 
is multiplied by the energy e^. 

4;SK 

Thus, our result should be compared to -37- ■ The critical 
compactness as a function of e-y for two values of the mag- 
netic field is shown in Fig. [TJ The two crosses mark the 
values of SK, which completely agree with ours in the limit 
of catastrophic losses. 

We note also that eq. @ , which was obtained by making 
robust calculations, also agrees with the eq. (OS!) in the limit 
of catastrophic losses. 



2.2.2. Step function approximation for the cr^^ cross section 

In this paragraph we analytically find the critical luminosity 
compactness after adopting a more realistic expression for 
the annihilation cross section: 



of the parameters. Thus, the scries of approximations con- 
verges and rig can be found in closed form: 



„ 00 

<(7) - ^VA"- = 



c 



7 



2 Z-^ 



(31) 



The luminosity in soft and hard photons can then be found 

(32) 



, ^ 4/b Ce^ 
T 3 2(1 -A) 



, _ 2^bC4 |a| 
T " 3 1-A' 

The critical 7-ray compactness is given by 

1 



/inj _ 1 

7,cr n 



3/2 



(33) 



(34) 



cttcto- 



e(a;e^ - 2) 



(26) 



Figure [D shows /"■j.j. as a function of for two values of 
the magnetic field. When the threshold energy for absorp- 
tion 2/e-y is much lower than the maximum energy of the 
soft synchrotron photons 6e^/4, the critical compactness as 
a function of Cy behaves as for the (5-function cross sec- 
tion. However, when the two energies become comparable, 
ll^^^j. abruptly increases, because the number of soft photons 
serving as targets for the absorption significantly decreases. 
This effect can be seen in Fig. [2] as a change in the curvature 
of the curves. 



3. Numerical approach 

To numerically investigate the properties of quenching, one 
needs to solve again the system of eqns (5)-(7) augmented 
to include more physical processes. As in the numerical 
code, there is no need to treat the time-evolution of soft 
photons and 7— rays through separate equations, the sys- 
tem can be written: 



where the normalization used is ctq = 4/3. The initial equa- 
tions ([S])-© keep the same form. Only the operator of pho- 



dne{'j,t) rio 



dt 



q: 



77 



syn 



(35) 
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Fig. 2. Log- log plot of the critical luminosity compactness 
as a function of the injected 7-ray photon energy e-y, 
wiien the 0-function approximation for the cross section 
cr-y-y was used for two values of the magnetic field strength: 
B = 40 G (soUd Unc) and B = 4 G (dashed hne). 



and 



dt 



7, esc 



V 

-'ICS 



(36) 



where and n-^ are the differential electron and photon 
number density, respectively, normalized as in §2. Here 
we considered the following processes: (i) Photon-photon 
pair production, which acts as a source term for elec- 
trons (Q°^) and a sink term for photons (>C^^); (h) syn- 
chrotron radiation, which acts as a loss term for elec- 
trons (i^syn) and a source term for photons (QJyn); (iii) 
synchrotron self-absorption, which acts as a loss term for 
photons (>C2sa) and (iv) inverse Compton scattering, which 
acts as a loss term for electrons (/Cf^s) and a source term 
for photons (Qj^g)- In addition to the above, we assume 
that 7— rays are injected into the source through the term 
(Qjnj)- The functional forms of the various rates have 



been presen ted elsewhere - Mastichiadis fc KirkI (jl997f) and 



iPetropoulou fc Mastichiadis! (|2009() . The photons are as 
sumed to escape the source in one crossing time, therefore 

^■7, esc — ^■ 

The parameters of the problem are completely specified 
once the values of radius R and of the magnetic field B are 



set and the injection rate Q^nj 



is specified. Because Q/nj is 



readily related to the injected luminosity L™^ through the 
relation 



= {m^cyV dxx Q 



(37) 



where V is the volume of the source, the injection rate is 
specified once the injection compactness £"^'' and the func- 
tional dependence of Q^^^ on e-y is set. The electron physical 
escape timescale from the source te,esc is another free pa- 
rameter which, however, is not important in our case. Thus, 



we will fix it at value 



R/c. The final settings 



Fig. 3. Critical 7-ray compactness as a function of when 
(i) a^j and synchrotron emissivity are approximated by 
a Q- and a d- function respectively and (ii) when for both 
quantities the full expressions are used. For the first case the 
numerical result (points) and the analytical one given by 
eq. (solid line) are shown. The full problem can only be 
treated numerically and the corresponding result is shown 
with a dashed line. The size of the source is i? = 3 x 10^^ 
cm and the magnetic field strength B = 3.57 G. 



are the initial conditions for the electron and photon num- 
ber densities. Because we are investigating the spontaneous 
growth of pairs and synchrotron photons, we assume that 
at t = there are no electrons or photons in the source, so 
we set ric(7, 0) = nj{x, 0) = 0. 

3.1. Monoenergetic j~ ray injection 

As a first example we study the case where the 7— rays 
injected are monoenergetic at an energy e^. The expected 
trivial solution of the kinetic equations is that because there 
are no soft photons in the source, the injected 7— rays will 
freely escape and their spectrum upon escape will simply 
be their production one. Electrons and photons at energies 
different from e-y will remain at their initial values, i.e., 
no(7,i) = n^{x,t) = 0. 

Before numerically examining the stability of the above 
system using the full expressions for the photon-photon an- 
nihilation cross section and the emissivities, we perform a 
test by comparing the analytical result given by eq. (1341) 
with the respective numerical. The two solutions completely 
agree as shown in Fig. [31 In the same figure we also plot 
the critical 7-ray compactness found numerically when all 
quantities used were given by their full expressions. The 
qualitative behaviour of the full problem solution is quite 
well reproduced by the analytical one obtained while using 
the 0- and 6- approximations for the cross section and syn- 
chrotron emissivity respectively. This allows us to use the 
analytical expression given by eq. (j34p for the full prob- 
lem as well, when investigating the properties of quenching 
at least qualitatively. From this figure it becomes apparent 
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that 7-ray quenching can affect all photons with energies 
that satisfy the feedback condition. Starting around this 
energy and moving progressively to higher ones, we find 
that the critical compactness first sharply decreases until 
it reaches a minimum and then starts increasing again, as 
■ ■ 1/2 

ll^l.^. oc e-y' asymptotically. The exact value of the critical 
compactness also depends on the strength of the magnetic 
field. For a higher B-field the required minimum energy 
ej decreases, which can be easily deduced by eq.([2). Note 
also that because of the marginal relation (p4)) , the required 
critical injected compactness becomes higher. Thus, by in- 
creasing B, the corresponding curve is shifted towards lower 
7— ray energies and higher critical compactnesses, keeping 
its shape (see e.g. Fig. 

Figure H] depicts the results of a series of runs where 
we kept constant and varied For sufficiently low 
values of the latter parameter, all 7— rays indeed escape 
from the system without any impedance. If we call the 
compactness of the outgoing radiation at energy e^, then 
obviously, in this case, = i'^ = t^^, where is the 
compactness of the outgoing luminosity. However, above 
a critical value of £]^\ which we will call i^^^c^., the loop 
described above starts operating, transferring the 7— ray 
luminosity to softer radiation. Thus, for i]^^ > f™!,^, we get 
= ll^ + e'^ ^ £1^+ e^^„ because £1^ saturates - by £1^ 
we denote the compactness of the reprocessed luminosity, 
which appears at lower energies. This can be seen in Fig. 21 
where the dashed line curve represents the dotted line £1^ 

and the full line which is always equal to £™j because 
there are no sinks of photons present. For the particular 
example shown in Fig. 01 the values of the other parameters 
used are, i? = 3 x 10^^ cm, B = 40 G and = 2.3 x 10"^. 
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-4 -3.5 -3 -2.5 -2 -1.5 

lOgly 

Fig. 4. Plot of the compactness of emerging hard 
(dashed), soft £^ (dotted), and total radiated luminosity 
(full line) as a function of the injected compactness ('J^K 
The 7— rays are assumed to be monoenergetic, injected at 
energy = 2.3 x 10^. The size of the source is i? = 3 x 10^^ 
cm, while the magnetic field strength is B = 40 G. The 
emerging soft radiation is a result of the self-quenching of 
the system. 



-2 - 



-3 - 




-12 -10 -8 -6 -4 -2 2 4 6 

lOgEy 

Fig. 5. Multiwavelength spectrum for a choice of £™'' such 
that £1^ = £^. The other parameters are the same as in 
Fig. m Full line corresponds to the case where all rele- 
vant radiative processes were taken into account, while the 
dashed line corresponds to the case where only the two 
'core' processes were considered, i.e., synchrotron and 77 
interactions. 

Figure [5] shows the multiwavelength (MW) spectrum 
for the particular value of £™\ which results in equipar- 

tition between soft and hard luminosities, i.e., £1^ = £^. 
While the hard luminosity is simply the escaping monoen- 
ergetic radiation at e-y, the soft spectrum mainly consists of 
the synchrotron component produced from the secondary 
pairs, which were injected at energies 7p = 7e = £7/2. The 
peak of this component is at epj^ — 67^, while it has a spec- 
tral index of —1/2 (or -1-1/2 in vF^ units) as a result of 
the electron cooling to a 7^^ distribution below the energy 
of injection. There is also a synchrotron self-Compton (ssc) 
component, which is greatly suppressed by Klein-Nishina 
effects. As can be seen from the figure, the corrections in- 
troduced from the inclusion of inverse Compton scattering 
and synchrotron self-absorption are negligible, at least for 
our set of parameters. 

Figure El shows the time evolution of the different lu- 
minosity compactnesses for the example case presented in 
Fig. |S1 Although the dynamical, non-linear system of equa- 
tions describing the system (see for example system of eqns. 
(IT5I) ') is of third order, which in principle can lead to unsta- 
ble solutions or/and 'chaotic' behaviour (SK), the system 
quickly reaches a stable state after a few crossing times. 
Figure ini shows that although initially absent, a soft pho- 
ton population builds up, which leads to a decrease of V^, 
because 7-rays are being absorbed by the soft photons. 

3.2. Power-law injection 

As a next step we investigate the effects of quenching for a 
power-law injection of 7— rays, i.e., when the injection rate 
can be written as 

— Qo£-y"0(£7 ~ £7111111 )0(e7max ~ £7)7 (38) 
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Fig. 6. Time evolution of the soft (dotted line), hard 
(dashed line) and total (solid line) compactness for the 
same set of parameters as in Fig. [5] 

where Qo is a normalization constant depending on the 
value of a, the low- and high- energy cutoffs of the distri- 
bution. 

Before proceeding to more quantitative calculations, we 
remark that obviously, for a given B the feedback criterion 
for quenching to start will be that of eq. 1^ where should 
be replaced by e-ymax- Assuming that this condition holds 
and that the compactness of the 7— rays is above the critical 
one, the production of synchrotron photons with energy 
a^max = ^f^max/^ is guaranteed. This means that 7— rays 
with energy > e^min.eft should interact in 77 interactions, 
where 

^7inin,cff — 2/ Xniax — ^/^^7inax' ('^^) 

Thus, depending on the relation between e-ymin.cff and e^min 
we can distinguish two cases 

(i) £7111111 < €7111111, eft: In this case quenching will affect 
7-rays from e-^max down to e7niin,eff- For < e^rnin^eS the 
spectrum will be unaffected, i.e., it will be the same as the 
one at injection. 

(ii) e^min,eff < £7111111: In this case quenching will affect 
the whole 7— ray spectrum. 

The above two cases are exemphfied in Figures [7] and ID 
In both figures we assumed that 7— rays with a power-law 
form of slope a = 2 are injected in a source of radius R = 
3 X 10"'^^ cm and magnetic field strength i? = 40 G. The only 
difference between the two cases is in the lower and upper 
limits, which are e^min = 23 and e7max = 2.3 x 10^ for the 
case shown in Fig. [7] and e-^min = 230 and e-ymax = 2.3 x 10^ 
for the case shown in Fig. [8l 

In the case of Fig. [7]we have £7111111 < £7111111, off = 1.7x 10'' 
and, as can be seen from the figure, the power-law breaks 
above £701111, eS because of quenching - indeed, we find that 
the value of £7111111, off is in reality lower by a factor of 2 than 
the one given in eq. ([39)) . because in the numerical calcu- 
lations shown here we use the full expressions for the syn- 
chrotron emissivity and the cross-section for 77 and not the 



5— function approximations implied in the derivation of the 
above relation. It is interesting to note that as the 7— ray 
luminosity increases and thus the quenching becomes more 
effective, the break between the unabsorbed and absorbed 
part of the power-law becomes more abrupt. 

In the case of Fig. [51 on the other hand, the choice of 
the parameters were such as to make £701111, off comparable 
to £7min- Thus, according to the qualitative analysis pre- 
sented above, quenching should affect the whole 7— ray dis- 
tribution. Indeed, as can be deduced from the figure, as the 
luminosity increases, absorption is moving progressively to 
lower 7— ray energies, until it affects the whole distribution. 
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Fig. 7. MW spectra for a power-law injection of 7— rays 
of the form ^ Qo^^^ for E^min < £7 < £7max with 
£7min = 23 and £7max — 2.3 X lO"' for various values of the 
normalization Qo- This starts from Qo = 4 x 10"'^ and then 
exceeds its previous value by a factor of 2. The dashed line 
curves represent the runs without quenching. The source 
parameters are _B = 40 G and i? = 3 x 10^^ cm. 

We also investigated the relation between the slope of 
the injected (ainj) and the 'quenched' (oquo) 7-i'ay spec- 
trum. For that reason, we ran the code for different ainj, 
choosing each time suitable values of the source parameters 
to enable quenching. Figure [S] shows the result of the runs 
(solid line) and the slope ainj + 1 (dashed line) for compar- 
ison reasons. We find that aquo fx Q^in^^- Thus, quenching 
becomes more pronounced as the injected spectrum gets 
steeper. 

4. Application to the one-zone homogeneous 
models: The 3C 279 case 

The above has a straightforward application to the so-called 
one-zone homogeneous models, which are customarily used 
to fit the MW spectr a of Active Galactic Nuclei (AGN) 
(|Bottcher et all 120091: lAleksic et all l201lt) . According to 
these models, the 'source' contains electrons and/or pro- 
tons, which radiate through various processes and form a 
spectrum that can then be fitted to the observed data by 
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Fig. 8. Same as in Fig.[7]with e^min ~ 2.3 x 10^ and e^max 
2.3 X 10^ with Qo starting from Qq = 2.5 x 10"''. 
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Fig. 9. Slope of the quenched 7-ray spectrum as a function 
of the injected spectrum's slope (solid line). To facilitate a 
comparison the line ainj + 1 is also shown (dashed line). 



varying certain key paramet e rs - f or a detailed descrip- 
tion s ee iMastichiadis fc KirS (|1997D and iKonopelko et al.l 
(|2n03D . 

Alternatively, one can assume that the 7— rays are pro- 
duced by some mechanism and investigate the allowed pa- 
rameter space that allows the escape of these 7-rays. To 
perform this analysis, one does not require the existence 
of other softer photons that act as targets, but one simply 
has to look for the threshold conditions for photon quench- 
ing. In other words, the relevant question one has to ask is: 
Assume a spherical source of radius R containing a mag- 
netic Held B and producing a 7— ray spectrum that is mea- 
sured at Earth with a flux F^{e^). If the source is moving 
with a Doppler factor S with respect to us, what is the pa- 
rameter space of R, B, and S that aUows the 7— rays to 



escape? As we will show, this can give very robust lower 
limits for the Doppler factor 6. 

The luminous quasar 3C 279 is a good example. A re- 
cent comprehensive re view of observations can be found in 
iBottcher et"all (|2009D . Here we will mainly focus on the 
2006 campaign, which disco vered the source at VHE 7- 
rays, showing a high TeV flux (jMAGIC Collaboration et al.l 
[2008,), while the X-rays were at a much lower level. 

To apply the ideas discussed in the previous section, we 
follow the procedure below: 

We assume a spherical source of radius R moving with 
Doppler factor d with respect to us and containing a mag- 
netic field of strength B. We assume that 7— rays are pro- 
duced by some unspecified mechanism inside the source 
and, upon escape, they form the observed 7— ray spectrum. 
If f™-" is the compactness of the 7— rays (c.f. eq. [T]), we can 
relate it to the observed integrated flux F'"* by the relation 



S^R 



(40) 



' for i/"*?^ < 
7 v"^-"" = 



where Dl = 3.08 Gpc is the luminosity distance to the 
source. Here we used a cosmology with fi™ = 0.3, f^A = 0.7 
and Hq = 70 km s"^ Mpc"^ 

To model the TeV emission, we performed a fit to 
the MAGIC data, cor rected for intergalactic 77 absorption 
(IBottcher et al.ll2009D . of the form vF^ = Av' 
V < <ax where A = 1.3 x lO^^ JyHz, a = 
2 X 10^5 Hz and i/jj^ = lO^^ Hz. Then we injected 7-rays 
in the code of the form — Qoe^" for e^min < < e-ymax 
where the normalization is set by condition (jSTj with the 
integration taking place between 



e7min = hiy^^l/ {SmcC^) and 



hK^JiSm^c^). 
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Fig. 10. MW spectra of 3C 279 for R = lO^'^ cm, B = 40 G 
and S = 10°-8 (dotted line), lO^ o (dashed line), lO^-^ (dot- 
dashed line) and lO"'^'^ (solid line). In all cases the 7-ray 
injection would have produced the solid line in the absence 
of quenching. Black squares and bowtie represent the ob- 
servational data from February 2006. 
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Fig. 11. MW spectra of 3C 279 for R = lO^^ cm, S = lO^ ^ 
and B =.16 G (solid line), 1.6 G (dashed line) and 16 G 
(dotted line). In all cases the 7-ray injection would have 
produced the solid line in the absence of quenching. Black 
squares and bowtie represent the observational data from 
February 2006. 



Fig. [TO] depicts the resulting MW spectra for R — 10^^ 
cm, B = 40 G and J varying between 10°-^ to 10^ '' with 
increments of 0.2 in the exponent. Only the run with the 
highest value of 8 can fit the TeV data. The other three, 
even though they have parameters that could also fit the 
data, they do not because of quenching. It is important to 
repeat at this point that the 7— ray luminosity produced in 
these cases is the maximum possible that the system can 
radiate for the particular choice of i?, B, and 8. Any effort 
to increase the 7— ray flux by e.g. increasing t^'^ would only 
result in an increase of the reprocessed radiation. 

Fig. [TT] depicts a series of similar runs with R = 10^^ 
cm, S = 101-2 and B = .16 G (solid line), B = 1.6 G 
(dashed line) and i? = 16 G (dotted line). While for low 
B values there is no quenching, this begins to appear at 
higher values of the magnetic field. Note, however, that the 
value of i? = 1.6 G produces an acceptable fit, because 
the soft reprocessed radiation does not violate any observ- 
able limit (most notably the X-rays) and it can produce a 
marginally acceptable fit to TeV 7— rays. Based on this, one 
can define a limiting value of the magnetic field i3q,mx(-R, S) 
with the property that if B > i3q,mx, the quenching does 
not allow the internally produced 7— ray luminosity to es- 
cape. Therefore -Bq.mx is the highest allowed value that the 
i?— field can take for the observed flux of TeV 7— rays. 

Figure [12] depicts the locus -Bq,mx as a function of the 
size of the source R for various Doppler factors 6. Each line 
separates the parameter space into two regimes: Because 
of the quenching, a choice of R and B above the locus 
produces a far too low 7— ray flux that cannot fit the TeV 
observations. For example, if one assumes that S = 10 and 
R — 2.5 X 10"'^^ cm - so that it corresponds to a variability 
timescale of 1 day, the magnetic field cannot exceed 0.3 
G. However, owing to the steep dependence of the flux on 



6 (c.f. Eqn 40), the above condition is greatly relaxed for 
higher values of 6. 

We should also mention that Bq^mx not only depends on 
5 and R, but also on the upper cutoff e^max of the injected 
power-law in 7— rays. For the calculations performed here, 
we chose e^niax such as to match the highest observed fre- 
quency in the data set, i.e., 10^^ Hz. If one assumes that the 
internally produced 7— rays can extend to even higher ener- 
gies, then the value of i?q,mx is affected as well. Accordingly 
we found that if we were to repeat the same calculations 
with e-yniax ten times higher than the value we used, then 
quenching would affect the flux more and the lines of Fig. 12 
would have to shift upwards by a factor of 2. Therefore, 
Fig. 12 can be taken as a conservative lower limit on i?q,mx- 
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Fig. 12. Parameter space of the allowed values of the mag- 
netic fleld B as a function of the size of the source R for 
Doppler factors of (5 =10 (dotted line), 15.8 (dashed line), 
and 25 (solid line). The space above a line for a certain 8 
corresponds to combinations of R and B that cannot suc- 
cessfully flt the TeV observations of 3C 279. 



5. Summary and discussion 

We have examined the effects of photon quenching on com- 
pact 7— ray sources. Automatic photon quenching, first pro- 
posed by SK, is a non-linear network of two processes 
that are very common in high-energy astrophysics, namely 
photon-photon annihilation and lepton synchrotron radia- 
tion. We verified and expanded both analytically and nu- 
merically the results of SK and we showed that this net- 
work, when applied to compact 7— ray sources, can give 
limits on the Doppler factor, which depend on the source 
radius and on the magnetic field strength only. This is a 
theoretically derived limit and as such does not depend on 
the soft photon populations that might be present inside 
the source. 

It is important to note that photon quenching can be 
investigated only when the radiation transfer problem is 
solved self-consistently because it is a non-linear process 
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that involves redistribution of energy from the 7— rays to 
the lower energy parts of the photon spectrum. For this we 
used a system of kinetic equations that treats in the sim- 
plest case 7— ray annihilation and electron synchrotron ra- 
diation. Using suitable (5— function approximations for the 
photon-photon annihilation cross section and for the pho- 
ton synchrotron emissivity, we were able to verify the re- 
sults of SK in the limit where synchrotron losses can be 
considered as catastrophic (see Fig. [T|). On the other hand, 
using a more realistic step function for the cross section, 
we were able to find a very good agreement with numeri- 
cal results, which utilized the full photon-photon annihila- 
tion cross section. Furthermore, we extended our numerical 
study in the case where 7— rays are injected in the form of 
a power-law. We found that quenching can be responsible 
for wide spectral breaks, which become more pronounced 
as the injected 7— rays become steeper (see Fig. [2]). 

Quenching has some interesting consequences when ap- 
plied to existing 7— ray sources. For example, for the 2006 
MAGIC observations of the quasar 3C279, modelling with 
combinations of low values of the Doppler factor S and high 
values of the magnetic field strength B lead to strong ab- 
sorption of the injected 7— rays. The only way to avoid 
quenching is either to assume typical values of 6 (i.e., 
around 10) and low values of B or, alternatively, high values 
for both S and B - see Fig. [T^l 

We restricted our analysis to cases where 7— rays are in- 
jected into the source without being too specific in assigning 
a particular radiation mechanism for producing them. We 
preferred this approach because our aim was to calculate 
the effects of quenching using simple functional forms for 
the 7— rays and not to model the source in detail. For ex- 
ample, we could have used proton synchrotron radiation to 
fit the 3C279 TeV observations. However, this would have 
involved extra parameters for the radiating relativistic pro- 
tons and would have complicated the analysis beyond the 
level at which we would like to present it here. 

In our analytical treatment we ignored inverse Compton 
scattering, because its inclusion would have complicated the 
analysis. On the other hand, it was taken into account in all 
our numerical calculations, although it is greatly suppressed 
by Klein-Nishina effects (see Fig. 5). This process can still 
make an impact in cases where the soft photon compactness 
greatly exceeds the magnetic one (c.f. Eqn. 14). However, 
this has to be seen in conjunction with the overall fitting 
of the MW spectrum of the source and we avoided present- 
ing it here, although we took it into account in §4 when 
performing the analysis on quasar 3C279. 

Concluding, quenching has many far-reaching implica- 
tions for the modelling of compact 7— ray sources such as 
Active Galactic Nuclei and Gamma Ray Bursts. Because 
quenching operates in an autoregulatory manner by cou- 
pling electrons to photons, its effects on the photon spectra 
cannot be manifested if the modelling is done by assum- 
ing an ad-hoc radiating electron and/or proton distribu- 
tion. This could lead to parameter choices that lie within 
the quenching 'forbidden' parameter space, making them 
in essence invalid. 

On the other hand, if taken fully into account, quench- 
ing can give robust limits for important source parameters 
such as its size, magnetic field, and the Doppler factor for 
a given 7— ray flux. In the case of 3C279 for example, we 
found that to avoid quenching, the source should have a 
large Doppler factor {S > 10) whenever high B values are 



adopted. Similar constraints could be found for other AGNs 
with spectra extending to even higher energies (~ 1 — 10 
TeV). In these cases the conditions for quenching are more 
stringent, because they critically depend on the maximum 
energy of the 7-ray photons and are certainly worth inves- 
tigating further. 
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